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Three-neutron resonances are investigated using realistic nucleon-nucleon interaction models. The 
resonance pole trajectories are explored by first adding an additional interaction to artificially bind 
the three-neutron system and then gradually removing it. The pole positions for the three-neutron 
states up to J— 5/2 are localized in the third energy quadrant - Im(E)^ 0, Re(E)^ - well before 
the additional interaction is removed. Our study shows that realistic nucleon-nucleon interaction 
jy^ i models exclude any possible experimental signature of three-neutron resonances. 
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I. INTRODUCTION 



The possible existence of pure neutron nuclei is a long standing ambiguity in nuclear physics. The neutron-neutron 
(nn) scattering length is negative and rather large a nn ~ —18.59 ± 0.40 fm 0, indicating that this system is almost 
' bound in 1 5'o state. This value is actually the signature of a virtual bound state just wlOO kcV above the threshold. 
^T) , It is thus expected that by adding a few neutrons, one can end up with a bound multineutron state, as it happens in 
some fermionic systems, like 3 He atomic clusters Q- This is the reason for the recurrent turmoils appearing in the 
nuclear physics society 0, 01- However, the weakness of nuclear interaction in higher partial waves (namely P and 
D), in comparison with their centrifugal barriers, excludes the possibility of binding 'virtual' dineutrons together |5j. 
Recently, the non-existence of small bound multineutron clusters has been settled out theoretically 0, IE 01 @j- 
Nevertheless, the existence of resonant states of such nuclei having observable effects had not yet been fully eliminated. 

Indeed, in spite of numerous experimental and theoretical studies that exploit different reactions and methods, the 
situation concerning few-neutron resonances is not firmly established. One does not have clear ideas even for the 
' simplest case: the three-neutron compound. A nice summary on the three-neutron system status up to 1987 can be 
£ , found in @. A few more recent experimental studies have provided only controversial results. In |10I. the analysis of 
Ch ■ 3 He(ir~ , 7r + )3n process yielded no evidence of a three-neutron resonant state. The explanation [llj of double charge 
exchange differential cross sections in 3 He in term of a broad (E = 2 — Qi MeV) three-neutron resonance were recently 
criticized by a more thorough experimental study 0] ■ Nevertheless the latter study further suggested the existence 
of a wide 3 n resonance at even larger energies (E sa 20 — 20z MeV). Furthermore, very recent experimental results on 
£3 , 8 He(d, 6 Li)4n reaction have shown some discrepancies with what could be expected from phase space calculations, 
suggesting the existence of a resonant tetraneutron [l3[ . 

There exists several theoretical efforts to find 3 n and 4 n resonances. A variational study based on complex-scaling 
and simplified nucleon-nucleon (NN) interaction was carried out in [l4T | with the prediction of a 3 n resonance at 
E = 14 - 13i MeV for the J* = 3/2+ state. On the other hand, no 3 n nor 4 n resonances were found by Sofianos 
ct al. using MT I-III potential; only the existence of some broad subthreshold (Er < 0) states was pointed out. 
Realistic NN forces can provide however different conclusions. These models contain interaction in P- and higher 
partial waves, which - due to the Pauli principle - are a crucial ingredient in binding fermionic systems. The only study 
performed using realistic potentials was carried by Glockle and Witala [l^. These authors were not able to find any 
real three-neutron resonances. However due to some numerical instabilities, the full treatment of 3 n system was not 
accomplished and the conclusions were drawn based on the phenomenological Gogny interaction ^8|. Reference 0] 
is probably the most complete study of three neutron system. In this work, the trajectories for 3 n states with J ^3/2 
have been traced by artificially enhancing a rank-2 separable nn interaction. Our work is devoted to complement the 
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studies of references 0,EJ- We have explored all 3 n states up to J=5/2, fully relying on realistic AW interactions. 
Similarly to Ref. |l7j we have used two different methods, namely complex-scaling (CS) and analytical continuation 
in the coupling constant (ACCC), to trace the resonance trajectories. This allows us to check the reliability of our 
results, as well as to judge on the pertinence of the two methods used. 

[Calculations presented in this paper use ^ = 41.44 MeV-fm 2 as an imput for the mass of the neutrons.] 



II. THEORETICAL BACKGROUND 



Resonance eigenfunctions correspond to complex energy solutions of the Schrodingcr equation: 

z 

Hty res — Eres^ res j E res — Eres ~^^res- (1) 

Since physical resonance have positive energy real parts, e res > 0, the corresponding eigenfunctions are not square 
intcgrable. Nevertheless, by applying them a similarity transformation they can by mapped onto normalizablc states. 
That is: 

(SHS^ (,S* res ) = E res (S^res) (2) 

with S^res — > as r — > oo. 

Functions (s^ re .s\ are in Hilbert space, although $ res are not. The complex-scaling method is defined by means 
of the similarity operator 0, |2(j : 

S = e l6r £, (3) 

such that any analytical function f(r) is transformed according to: 

Sf(r) = f{re w ). (4) 

For a broad class of potentials, complex-scaling operation does not affect the bound and resonant state spectra of the 
Hamiltonian H, provided ^ 8 < |. However the continuous spectra of H will be rotated by an angle 26. Resonance 

eigenfunctions ^S^resJ of the scaled Hamiltonian become square integrable if 26 > |arg(£' res )| and therefore the 
standard bound state techniques can be applied to determine the corresponding eigenvalues. 

In order to solve the three-body problem we use Faddeev equations in configuration space |2l| . first derived by 
Noyes [2^ ■ Though this formalism was initially developed to investigate the three-body continuum, it turns to be 
very useful to treat bound-state problems as well. For three identical particles, the Faddeev-Noyes equations read: 

(e-H - %) ih 3 .k = V lJ (P+ + P-)tp ijt k. (5) 

Hq is the three-particle kinetic energy operator, Vij the two-body force, V'y'jfc the Faddeev component and P + ,P~ 
denotes cyclic particle permutation operators. The properly symmetrized three-body wave function is \& = (1 + P + + 
P~)ipij,k- In order to simplify the kinetic energy operator and to separate internal and center of mass degrees of 
freedom we use Jacobi coordinates: Xij = Tj — T\ and y^j = [r^ — i (j\ + r^)] . 

Complex-scaling Faddeev equations causes no difficulties: one has simply to scale all the Jacobi vectors with the 
same exponential factor: 

x^ — > x^e 10 and —> y^e 10 . (6) 



Such transformation affects only the hypcrradius p = J x 2 j + y 2 ^ , whereas neither the expressions of permutation 
operators - P + ,P~ - nor the angular dependence of Faddeev equations are affected. Using Jacobi coordinates, the 
kinetic energy operator can be expressed as a six-dimensional Laplacian Hq = — ^-A x with \ = (xij, y~ij) and equation 
© transforms into: 

[E + e- l28 A x - (xijj 6 )} $ij,k(xy,Vv) = Vij (x^e 10 ) (P+ +P"fe(^,^) (7) 
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As in our previous works [flU^, equation (J7J) results into a system of integrodifferential equations by first expanding 
and then projecting the Faddeev components into a - partial wave - basis of spin, isospin and angular momentum: 

T i — > — >\ S^^k,l x lS X ijiVio) r a T — An ,_, 

Vij.k [Xij , Vij ) = 2^ l e ®yi,iy {Xij , yij )\ JT . (8) 

STL Xi i y » 

l x ly 

Next, the radial dependence of amplitudes ip^Tj, ( x ij>Vij) 1S developed in a basis of cubic Hermite splines |24j . 
Such a procedure applied to the complex scaled Faddeev equation J7J results into a generalized eigenvalue algebraic 
problem: 

A X = E res B (9) 

where A and B are known complex matrices, E res and X are respectively the complex eigenvalue and eigenvector to 
be determined. 

The partial wave expansion JHJl is pushed up to amplitudes with intermediate angular momenta max(j x , j y ) <4.5, 
which guarantees at least three-digit accuracy in the results. This requires solving linear system 10 of a relatively 
large (n eq ~ 10 5 — 10 6 ) size, which prevent us applying direct linear algebra methods. In order to avoid cumbersome 
matrix inversion, we use inverse-iteration techniques to search only individual eigenvalues and iterative methods to 
solve the linear systems. Technical details of the numerical methods in use can be found in ||. 

There is an apparently simpler procedure to depict resonance trajectories, namely the method of analytic continu- 
ation in the coupling constant (ACCC). This method, developed by V.I. Kukulin et al. [H HI [H HI 113 , is based 
on the intuitive argument that a resonance may be treated as an eigenstate which arises from a bound state when 
the intensity of the attractive part of the interaction decreases. The S-matrix pole of a resonant state is defined as 
the analytic continuation of a bound-state pole in the coupling constant of the attractive part of the Hamiltonian. 
One assumes that the Hamiltonian can be written as H(X) = Hq + XH at t, where H att is the attractive part of the 
unperturbed Hamiltonian H(X = 1). If for some value of A the Hamiltonian has a bound state, then by gradually 
decreasing A, its binding energy decreases and the state reaches the threshold at A = Ao - i.e. E(\ = Xq) = - to 
become a resonant or a virtual one. 

For physical Hamiltonians, the binding energy is assumed to be analytical in A. Moreover, it can be shown H3j 
that for a two-body system the square root of the binding energy - kg = y/—E - behaves near the threshold A = Ao 
as: 



By introducing the complex variable 



A-A for £ = 
ke ~< VX—J- Q for £^Q ( 10 ) 



A - A for I = 



VA- A for £^0 



one can consider kg as a function of variable x and analytically continue it from the bound state (A > Ao) to the 
resonance re gion (A < A ). Motivated by the functional form (|1L)|I near the threshold, we use for kg the Pade 
approximant |30(: 

kl ( X ) = k^\ x) = - a r x+a r x2 2 + --- a h N ' xN M r in) 

w w 1 + h ■ x + b 2 ■ x 2 + . . . b M ■ x M ~ l K ' 

In general, for systems with n > 2, kg is the relative momenta to the nearest desintegration threshold kg = 
^—{E n — Ei <n ). Contrary to the n=2 case, for n > 2 systems the angular momenta does not determine anymore if 
a bound state turns into a virtual or resonant one when A < Ao- Nevertheless, according to l|10f) this transition can 
be discriminated by studying the nearthrcshold behavior of ^/—(E n — Ei <n ). If (E n — Ei <n ) is linear in (A — Ao) at 
the origin, the bound state turns into a resonance; virtual state appears if (E n — E i<n ) is quadratic in (A — A ). 
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III. THREE-NEUTRON RESONANCES 



Complex scaling methods (CS) turns out to be a very powerful tool when treating resonance problems in atomic 
and molecular systems. However scaled nuclear potentials introduce numerical instabilities not encountered in atomic 
physics, which is dominated by Coulomb (or Coulomb derived) potentials. The long range part of the realistic 
nuclear potentials decreases exponentially with a Yukawa pion-tail. Using CS, exponentially decreasing function e~ fJ,r 
transforms into e -f J -l cose + lsme ] r . For a scaling parameter 9 > 0, the potential range increases as l/cos9, introducing 
sizeable oscillations and demanding larger and denser grids to describe the system. These difficulties can be partially 
avoided by extending 6 to complex values, as proposed by Glockle et al. [l6j . 
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FIG. 1: Reid 93 and AV18 l So nn potentials. 

A much more serious problem arises due to the hard-core short range part of the nuclear potentials (see Figs.^). 
Notably, the presence of the e~ c ^ terms in AV18 Hi and AV14 ^ models settles that for CS with 9 > 45° these 
potentials become divergent. Furthermore even for smaller 9 these two models have a strong oscillatory behavior 
which makes numerics unstable. This can be seen in Fig. where the transformed NN potentials for two different 
angles are plotted. Note that the number of oscillations as well as their amplitudes increase dramatically with the 
rotation angle. In this sense Reid 93 model (3(| has the best analytical properties with respect to CS, though it 
provides stable results only for the 9 < 35°. 

When treating two-body systems, the numerical complications described above can be avoided by using the so-called 
'smooth-exterior' complex scaling (SECS) 0. The problem due to the hard core of the nuclear potentials is overcome 
by performing scaling operation only in the outside region of the interaction. This method can be implemented by 
using the following transformation UM '■ 



with 



F(r) 



¥ 6 - 1] 



F(r) 



I [1 + e 2»7(r-fl)] [1 + e -2^(r-fl)] 

4^ ln [1 + e 2 ^ r + R )] [1 + e- 2r ^+ B ^} 



(12) 



However, despite the effici ency of this method when dealing with n=2 

lems t 
t 0.5. 



where ry and R are smoothing parameters. However, despite trie efficiency 

problem, it is not easy to implement it in the n=3 Faddeev equations |l6l |. The problems turned out to be crucial 
when calculating deeply lying resonances and we were limited to those with E re /Ei mg 
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FIG. 2: Analytical continuation of the 1 So Reid 93 and AV18 nn potential under the complex scaling transformation eq. (J1J 
with 8 = 20° and 8 = 45°. One can notice that AV18 potential results into much more oscillating structure than Reid 93. 



All the difficulties discussed above are not encountered in the ACCC method, which can be used in principle even 
to determine subthreshold (s res < 0) resonances. In calculating the input binding energies for determining the Padc 
approximant (|llf> . one can also scale the total nn potential with a scaling factor 7 - not necessarily only its attractive 
part - and consider 7 as the extrapolation parameter 7 = A. 
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FIG. 3: Comparison between ACCC and CS methods for 3 Po nn resonance trajectories with Reid93 potential. ACCC results 
with several Pade orders [N,M] and 7=6.1 to 1.0 are represented by solid lines. CS results are denoted by squares and correspond 
to 7 from 6.1 to 2.7 by steps of 0.1. Star points correspond to [5,5] Pade approximant used in ACCC. They are already very 
close to CS results: by adding a few terms in the approximant a perfect agreement is reached. 



In Figure |3 we compare the results of various order Pade approximant used in (|11|) with the ones obtained using 
CS method: dineutron resonance trajectory is depicted for Reid 93 3 Pq rm-waves. One can see that a nice agreement 
is reached between both techniques even when calculating broad resonances, near and beyond the saturation point 
(where real energy part starts decreasing), provided high order Pade approximant are used. This requirement can 
not always be met in numerical calculations with realistic interactions, since it implies a very accurate binding energy 
input. Another critical point of ACCC method, as discussed in [23, is that its efficiency highly depends on the 
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accuracy of Ao value. The determination of Ao is often difficult, since one has to deal with critically bound systems 
having very extended wave functions. Therefore one has to be rather prudent when applying this method and always 
check if Pade extrapolation converges. 

TABLE I: Critical enhancement factors 7 C required to bind dineutron in various states and for different 7V7V realistic interaction 
models in use. 





Nijm 


II Reid 93 AV14 AV18 


M L s ) 


1.088 


1.087 


1.063 1.080 


2 n( 3 P ) 


5.95 


5.95 


5.46 6.10 


2 n{ 3 PF 2 ) 


3.89 


4.00 


4.30 4.39 


2 n{ 1 D 2 ) 


9.28 


9.22 


9.54 10.20 



Before starting to analyze three-neutron system it could be useful to discuss the basic properties of dineutron and 
nn interaction in general. As mentioned above, dineutron is almost bound in the 1 Sq state: one should enhance the 
nuclear potential only by a factor -f c w 1.08 to make it bound (see Table QJ. Notice the very good agreement for the 
critical enhancement factors 7 C obtained by different local AW-interaction models. Only AV14 slightly deviates from 
the other models, due to its charge invariancc assumption: the potential being adjusted to reproduce neutron-proton 
(np) scattering data, ignores the fact that experimental So nn scattering length is smaller in magnitude than np 
one 0. The spherical symmetry of this state makes that when 7 — > 1, i.e. for the physical value of the potential, the 
bound state pole moves staying on the imaginary fc-axis and becomes a virtual state, not a resonance. The approximate 
position of this state can be evaluated from the nn scattering length by means of the relation E virt w ma2 . The 
results thus obtained together with the exactly calculated virtual state energies are given in Table ITT1 

TABLE II: Nuclear model predictions for nn scattering length (in fm) as well as corresponding virtual state (in MeV), evaluated 
from scattering length and calculated exactly. 





Nijm II Reid 93 AV14 AV18 


flnn( 1 5'o) 

h 2 /(mal n ) 

E v i r t{ So) 


-17.57 -17.55 -24.02 -18.50 
0.134 0.135 0.072 0.121 
0.1162 0.1165 0.0647 0.1055 



In fact, multineutron physics being in low energy regime, is dominated by the large nn scattering length value 
(a nn )- The wave function has only a small part in the interaction region (effective range ro << a nn ) and therefore 
marginally depends on the particular form that the nn 1 5o-potential can take, once ro and a nn are fixed [f|- On the 
other hand ro is controlled by the theory - the pion Compton- wavelength - and a nn is constrained by the experiment. 
These effective range theory arguments [3lJ show that one should not rely on the modifications of 1 £0 waves to favor 
the possible existence of bound or resonant multineutron states. 

P-waves of nn interaction are extremely weak and this turns to be a major reason why multineutrons are not 
bound 0. Neutron-neutron interaction in 3 Pi channel is even repulsive, whereas potentials in 3 PF2 and 3 Po channels 
should be multiplied by considerably large factors - 7 = [3.9 — 4.4] and [5.5 — 6.1] respectively (see TablcQJ - to force 
dineutron's binding. The centrifugal barrier moves these artificially bound states into resonances when 7 is slightly 
reduced from its critical value. 

Dineutron resonance trajectories for 3 Po and 3 PF2 states are displayed respectively in Figures 0] and [SJ They 
have been obtained with CS and ACCC methods and different NN interaction models. By using a high order Pade 
approximant and accurate 2 n binding energies input in the ACCC method, we observe a perfect agreement between 
these two techniques. Due to the limitations in the rotation angle 9, CS method was applied only up to 7 = 2.7. 
The corresponding resonance positions obtained with different NN models are explicitly indicated in the figures. One 
should note that the resonance trajectories in both states have similar shapes. When decreasing 7 from its critical 
value 7c, the resonance width starts departing slowly from zero and then increases linearly On the other hand, the 
real part of the resonance energy - e res - first grows linearly with 7 C — 7. Afterwards, this growing saturates and 
reach a maximum value from where it quickly decreases, vanishes and becomes negative. ACCC method allows to 
follow the resonance trajectories up to the physical value 7 = 1 at which the dineutron states are deeply subthreshold. 
However, the transition to the third energy quadrant happens at a value 7 C well above. Some resonance trajectory 
properties obtained using ACCC method are summarized in Table IIIII Note that in none of Figures Q] and [5] the 
resonance trajectories are plotted until the physical value 7 = 1. 

It is clear from results given in Tablc lrTTl that the existence of any observable P-wave dineutron resonance is excluded: 
only subthreshold poles with large widths persist, making such structures physically meaningless. It is worth noticing 
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recent few-nucleon scattering calculations indicating that a good description of 37V and 47V scattering observablcs, 
would require stronger 7V7V P- waves 133 . l33l | . However the modifications involved - less than 20% - let the dineutron 
resonances still subthreshold 



(MeV) 



nn 




FIG. 4: Dineutron 3 Po resonant trajectories in complex energy plane for AV14, Nijmll, Reid 93 and AV18 nn interactions. 
Solid dot symbols correspond to 7 values from 6.1 to 2.7 by steps of 0.1 obtained with CS method and Reid 93. Continuous 
lines represent ACCC results. ACCC and CS results superimpose up to 7=2.7 point, the limit of CS applicability. 



TABLE III: Enhancement factor values (7^) at which dineutron resonances become subthreshold (e re s = 0), and imaginary 
energy values Ei m9 (7c) at this point (MeV). Resonance energy E res for physical nn interaction (i.e. at 7 = 1) obtained using 
ACCC method. 





»P 


S PF 2 




Nijm II 


Reid 93 AV14 


AV18 


Nijm II 


Reid 93 AV14 


AV18 


7c 


2.27 


2.26 2.08 


2.24 


1.64 


1.71 1.46 


1.73 


Eim 9 (7c) 


-10.2 


-10.3 -10.6 


-10.2 


-45.6 


-36.9 -56.2 


-40.3 


E res (7 = 1) 


-14.l-17.2i 


-14.2-18. 5i -10.3-18. li 


-12.1-18.0j 


-20.5-64.8i 


-15.9-39.9j -17.9-80.li 


-34.l-45.4i 



One should remark the astonishing similarity for the P-wave dineutron resonance trajectories obtained with NN 
potientials having quite a different shapes (see Fig. 3 P curves for the three charge symmetry-breaking models 
considered superimpose, whereas for 3 PF2 they separate only when very large resonance energies are reached. En- 
hancement factors employed in tracing these curves are unphysically large and produce very broad resonances: 3 Po 
state moves into third energy quadrant at Ei mg ~ 10 MeV, while in 3 Pi<2 case this value goes beyond 35 MeV. 

Two neutron states with orbital angular momentum £=2, can be realized only in singlet state ( J £>g). This state is 
dominated by a sizeable centrifugal barrier and the critical enhancement factors to bind dineutron (7 C ) is considerably 
large (see Table QJ. The central effective potentials 



V eff (r) = V nn (r) + 



h 2 l{l + 1) 



in this and higher angular momentum nn partial waves are smoothly decreasing functions, without any dips, implying 
without further calculation that dineutron can not have I ^2 observable resonances in the 4th quadrant. 
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As mentioned in section [HI we cannot obtain numerically all the eigenvalue spectra of the 3n problem. Only a few 
specific eigenvalues of the linear system (0 can be extracted by applying iterative methods. When using the CS 
method, these techniques do not allow us a-priori to separate eigenvalues related to the resonances from the spurious 
ones related to the rotated continuum. In order to force the numerical procees to converge towards the resonance 
position we have to initialize it with a rather accurate guess value. This obliged us to follow the procedure described 
in [Ti| : first we bind three neutrons by artificially making tin interaction stronger, and then we gradually remove 
the additional interaction and follow the trajectory of this state. Note, that in bound state calculations one can use 
linear algebra methods determining the extreme eigenvalues of the spectra (e.g. Lanczos or Power-method), whereas 
resonance eigenvalue in the CS matrices is not anymore an extreme one. 

By enhancing the 1 So ran-potential there is no way to bind three neutrons without first binding dineutron. On the 
other hand, as quoted before, this wave is controlled both by theory and experiment and modifications of its form can 
not affect multineutron physics. Three-neutron can neither be bound if we keep the 1 So interaction unchanged, and 
multiply all nn P-waves with the same enhancement factor; by doing so, dineutron will be bound in 3 Pp2 channel 
before three-neutron was. In view of that, we have tried to enhance only one of P interactions, keeping the usual 
strengths for the other ones. The 3 Pi potential is purely repulsive and its enhancement can not give any positive 
effect. The Pq enhancement gives null result as well: dineutron is always bound before any 3 n states is formed. 
Only enhancing 3 PP2 channel we managed to bind trineutron without first binding dineutron and this happens in 
the J 71 " = | 3 n state alone. These properties turned to be general in the four realistic interactions (AV14, Reid 93, 
Nijm II and AV18) we have investigated. 

The critical 3 PP2 enhancement factors 7 C required to bind 3 n are summarized in Table IIVI They are so large 
that dineutron, although unbound, is already resonant; the corresponding resonance positions are summarized in 
the bottom line of this Table. One can remark, once again, the rather good agreement between the different model 
predictions. This fact, as well as the similarity in the dineutron case, predictions suggest that the different realistic 
local-interaction models would provide a good qualitative agreement in multineutron physics as well. Therefore, in 
further analysis of three-neutron resonances we will rely on a single interaction model. In this respect, Reid 93 model 
is the most suited, since it possess the best analytical properties and consequently provides the most stable numerical 
results for CS method. 

As mentioned, above only a | three-neutron state can be bound by enhancing single NN interaction, without 
first binding dineutron. In Figure we trace the 3 n resonance trajectory (full circles) for this state when the 3 Pi<2 
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TABLE IV: Critical enhancement factors 7 C required for 3 PF2 nn channel to bind J" = | three- neutron and corresponding 
J w = 2~ dineutron resonance positions (MeV). 





Nijm II 


Reid 93 


AV14 


AV18 


7c( d n) 


3.61 


3.74 


3.86 


3.98 


E( 2 n) MeV 


5.31-2.4H 


5.41-2. 52i 


5.20-2.49i 


4.83-2.3H 



enhancement factor changes from 3.7 to 2.8 with step of 0.05. Result were obtained using CS methods. Extending the 
calculations to smaller 7 values generated numerical instabilities, due to the necessity to scale Faddeev equations with 
ever increasing 8. It can be seen however, that this trajectory bends faster than the analogous one for the dineutron 
in 3 Pi<2 state, therefore indicating that it will finish in the third energy quadrant with Rc(E)<0. 

Three-neutron can also be bound in states | and ^ states by enhancing combined 3 PF2 and 3 Pi waves. However 
such a binding is a consequence of strongly resonant dineutrons in both mentioned waves. These resonances are very 
sensible to the small reductions of the enhancement factor and they quickly vanish leaving only the dineutron ones. 




FIG. 6: = 3/2~ three-neutron state resonance trajectory obtained when reducing the strength W of phenomenological 
Yukawa-type force (open circles for CS and solid line+star points for ACCC methods). Trajectory depicted by full circles 
represents one obtained using CS, when reducing enhancement factor 7 for 3 P2 — 3 F2 nn interaction . Trajectory depicted by 
full squares is dineutron resonance path in 3 P2 — 3 F2 channel, obtained by enhancing nn-interaction in these waves. Presented 
results are based on Reid 93 model. 



In order to systematically explore the possible existence of resonances in all the three- neutron states, we let un- 
changed the original AW-interaction and force the binding by means of a phenomenological attractive three-body 
force. We have assumed the latter to have an hyperradial Yukawa form: 

p_ 

V 3n = -W— , with p=yj 4 + y\. (13) 

with pq = 2 fm. In this way, dineutron physics is not affected. 

In table [V] we summarize the critical values Wo of the strength parameter W for which the three- neutron system 
is bound in different states. Corresponding resonance trajectories obtained by gradually reducing parameter W are 
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TABLE V: Critical strengths Wo in MeV-fm of the phenomenological Yukawa-type force of eq. (1131 required to bind three- 
neutron in various states. Parameter po of this force was fixed to 2 fm. W' are the values at which three-neutron resonances 
become subthreshold ones, whereas Btrit are such 3NF corresponding triton binding energies in MeV. 
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traced in Fig. [7] As in previous plots, CS results are presented by separate solid points, whereas ACCC ones are using 
continuous line and star points. One can see a very nice agreement between the two methods except for the J* = | + 
state, where a discrepancy between two methods manifests at large energy. This is probably an artifact of the very 
strong 3NF used. Such a force confines the three-neutron system inside a «1.4 fm box, a distance smaller than the 
3NF range itself, and starts to compete with the hard-core repulsive part of the nn interaction, making the ACCC 
method badly convergent for broad resonances. For the J n = | + state, due to the even larger W values, ACCC 
method has not been used. 



(MeV) 




FIG. 7: Three-neutron resonance trajectories obtained when varying the strength W of Yukawa-type 3NF. Results obtained 
using CS method are plotted with solid points. For the = 1/2 - state, W varies from 520 to 300 MeV-fm by steps of 20 
MeV-fm, denoted for shortness [520,300,20]. For ,F = 3/2", W is in the range [420,270, 10] and for ,F = 1/2+ in [310,210,10]. 
Other 3n states - depicted in the smaller figure - require considerably stronger 3NF to be bound and result in large resonance 
energies values. For J 71 = 5/2" W changed in [610,450,10], for ,F = 5/2+ W in [810,750,10] and then in [750,510,20], and for 
,r = 3/2+ in [1060,800,20] and [800,640,40]. ACCC method results are presented by solid lines, supported by star points. 



One can remark that resonance trajectories have similar shapes for all the three-neutron states, with the resonance 
poles tending to slip into the third quadrant (e res < 0) well before W equals 0, i.e. when the additional 3NF are fully 
removed and only ATV-interaction remains. In Table IVl are also given the values W' , obtained using ACCC method, 
at which 3 n resonances become subthreshold (e res = 0). These values are still rather large, strongly exceeding what 
could be expected for realistic 3NF. To illustrate how strongly such 3NF violate the nuclear properties, we give in the 
last row of the Table, the triton binding energies, obtained supposing that the same 3NF with strength W' acts on 
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it. These energies would be even larger if the range parameter po taken in our 3NF model had had smaller and more 
realistic values. 

The preceding results demonstrate that realistic ATV-interaction models exclude the existence of observable three- 
neutron resonances. In 0] 3 n, resonance in | + state was claimed at E = 14— I3i MeV for the non-realistic Minnesota 
potential. Our results using realistic nn interaction contradict however the existence of such a resonance. A very 
strong additional interaction is required to bind three-neutron in J v = | state. Removing this interaction, the 
imaginary part of the resonance grows very rapidly while its real part saturates rather early (starting with sal060 
MeV-fm, it reaches its maximal value at W «720 MeV-fm). Once this saturation point is reached, the resonance 
trajectory moves rapidly into the third quadrant. 

The results shown in Figure represent the 3 n resonance trajectories only partially, without following them to their 
final positions W = 0, i.e. when the additional interaction is completely removed. The reason is that these positions 
are very far from the bound state region, thus requiring many terms in Pade expansion to ensure an accurate ACCC 
result. In principle, one could always imagine that these trajectories turn around and return to the fourth quadrant 
with positive real parts. Although we have never encountered such a scenario in practical calculations, we are not 
aware of any rigorous mathematical proof forbiding it and it cannot consequently be excluded. 

Such a kind of evolution seems however very unlikely in a physical case of interest. The resonance trajectories 
should indeed exhibit a very sharp behavior once entered in the third quadrant whereas all our results indicate always 
a rather smooth variation. 

In this study we have deliberately omitted the use of realistic 3NF models. The reason is that such forces are not 
completely settled yet, specially for pure neutron matter. In addition, we should remark that the UIX [53 3NF acts 
repulsively for multineutron systems The more recent Illinois 3NF contains charge symmetry breaking (CSB) and 
considerably improves the underbinding problem of neutron rich nuclei present for AV18+UIX |3Sj . However, even 
strongly CSB realistic 3NF would be by order too weak to make three-neutron system resonant. The reason of weak 
3NF efficiency in multineutron physics is that such force requires configurations when all three neutrons are close to 
each other, whereas such structures arc strongly suppressed by Pauli principle. 

Finally, one can expect that the different ways used to artificially generate bound states could give rise to different 
resonance trajectories. Some resonance could thus have existed which are missed in our approach. To investigate 
such a possibility we have chosen a J 77 = | resonance, obtained by means of the phenomenological 3NF force <|13[) 
with W = 360 MeV-fm. Then we gradually reduce W to zero, increasing at the same time the enhancement factor 
for the 3 PF2 potential from 7 = 1 to 7 = 3.7. The resonance trajectory obtained this way is plotted (cross circles) in 
Figure |H1 together with the resonance curves obtained by additional 3NF only (open circles) and by enhancing only 
3 PF2 nn potential (full circles). Once 3NF was completely removed, the resonance pole joined the curve obtained with 
enhanced 3 PF2 channel. Note that the structure of bound state obtained with 3N force and by enhancing P- waves 
are quite different. 3NF requires very dense and spherical symmetric neutron wave functions. This is the reason why 
the | + state is more favorable than | (see Table 0) to bind three-neutron when using such an additional force. 

IV. CONCLUSION 

A systematic study of three-neutron resonances using realistic AW-interaction models was presented. The search 
of resonance positions was carried out by artificially enhancing the interaction between the neutrons in such a way 
that the three-neutron system becomes bound. Then, by gradually removing the additional interaction we followed 
the path of the resonance energy located in the third and fourth quadrants of the second Riemann sheet. 

Two different methods were successfully applied to trace resonance trajectories, namely the Complex Scaling (CS) 
and the Analytical Continuation in the Coupling Constant (ACCC). They provided results in very good agreement. 
Two alternative ways were also explored to enhance the interaction between neutrons: in one of them, the three- 
neutron system was bound by adding a phenomenological threc-nucleon force (3NF) and in the other one the binding 
was obtained by enhancing the interaction in some NN partial waves. 

All 3 n resonance trajectories, for states up to J = 5/2, were shown to move into the third energy quadrant 
(Re(E)< 0) becoming subthreshold resonances, well before the additional interaction is fully removed. Our results 
clearly demonstrate that the current realistic interactions exclude the possible existence of bound as well as resonant 
three-neutron states. To push the three-neutron resonance out of the subthreshold region, such models should already 
be strongly violated. These findings support the results of Ref ^} for simplified AW-interaction model. 

The possible existence of observable four-neutron (tetraneutron) resonance seems rather doubtful as well, since 
for such a system to be artificially bound one requires almost as large enhancement factors as in the three-neutron 
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FIG. 8: J w — 3/2 3n resonance trajectories obtained when reducing the strength W in the Yukawa-type 3NF (open circles) 
and enhancing 3 PF2 nn interaction (full circles). Crossed circles indicate the resonance path obtained when starting from 
(W = 360 MeV-fm,7 = 1) and going gradually to (W = 0,7 = 3.7). 

case Still, such a possibility can not be completely excluded, since the tetraneutron is favored due to presence of 
two almost bound dineutron pairs. Moreover, recent experiment on s He(d, Li)4n reaction shows an excess of low 
energy 6 Li nuclei with respect to what one could expect from a phase space analysis ^3|- The presence of events 
associated with a resonant tetraneutron state was suggested. The possible existence of such structures will be explored 
in a forthcoming work. 
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